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Abstract. We present a detailed study on the possibility of manipulating quantum 
information encoded in the "radial" modes of aiTays of trapped ions (i.e., in the ions' 
oscillations orthogonal to the trap's main axis). In such systems, because of the tightness 
of transverse confinement, the radial modes pertaining to different ions can be addressed 
individually. In the first part of the paper we show that, if local control of the radial 
trapping frequencies is available, any linear optical and squeezing operation on the locally 
defined modes - on single as well as on many modes - can be reproduced by manipulating 
the frequencies. Then, we proceed to describe schemes apt to generate unprecedented 
degrees of bipartite and multipartite continuous variable entanglement under realistic noisy 
working conditions, and even restricting only to a global control of the trapping frequencies. 
Furthermore, we consider the transmission of the quantum information encoded in the radial 
modes along the array of ions, and show it to be possible to a remarkable degree of accuracy, for 
both finite-dimensional and continuous variable quantum states. Finally, as an application, we 
show that the states which can be generated in this setting allow for the violation of multipartite 
non-locality tests, by feasible displaced parity measurements. Such a demonstration would 
be a first test of quantum non-locality for "massive" degrees of freedom (i.e., for degrees of 
freedom describing the motion of massive particles). 
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1. Prologue: the promise of alternative continuous variable degrees of freedom 

The last decade saw a boom in the development of experimental capabilities available for 
quantum information processing. The ability to manipulate the information of discrete 
variables encoded in polarisation, spin and internal atomic degrees of freedom has by now 
reached very high standards. On the other hand, the control and manipulation of continuous 
variable (CV) quantum information is still almost exclusive to light fields in quantum optical 
settings. Even though purely optical systems rely on well established tools and are the 
natural choice for communication tasks over long distances, they also suffer from significant 
drawbacks. Notably, the entanglement generation in such systems is strongly limited by the 
efficiency of parametric processes in nonlinear crystals; moreover, 'static' optical degrees of 
freedom - i.e. light resonating in cavities - are seriously affected by losses and decoherence 
over their typical dynamical time scales. 

Such limitations motivate the question of whether the full potential of infinite 
dimensional Hilbert spaces could be better harnessed by "massive" {i.e. , related to the position 
of a massive particle) CV degrees of freedom. Of course, to compete with the so far very 
successful quantum optical toolbox, such degrees of freedom would have to allow for a 
range of coherent manipulations at least as exhaustive as what quantum optics currently 
permits. Besides, to be considered advantageous over quantum optical systems, such degrees 
of freedom should allow for notable improvements in the generation of quantum entanglement 
and squeezed states under realistic working conditions. 

In the present paper, we argue that trapped ions meet both such requirements, and present 
an extensive study to substantiate this argument. In particular, we shall focus on the radial 
motion of trapped ions, that is on oscillations along a direction orthogonal with respect to 
the array of ions. These oscillations are described by continuous variable quantum degrees 
of freedom which we will refer to as radial modes lU [JJ. Radial modes have attracted 
considerable interest in the last few years, mostly in view of the fact that they allow for a 
tighter confinement (which also permits one to define the phonons locally). 

The paper is organised as follows. After having introduced the description of the physical 
system (Sec. |2]i, in Sec. [3] we demonstrate that any linear optical and squeezing operation 
can be obtained for radial modes of trapped ions by controlling the individual radial trapping 
frequencies, indicating that trapped ions can at least match the processing capabilities possible 
for hght modes. In Sec. |4]we show that, even restricting to cases where only global control 
of the trapping potentials is possible, such systems are actually apt to outperform optical 
modes in the generation of entanglement, both bipartite and multipartite. As applications we 
consider, in Sec. |5] the propagation of quantum information along the array of ions, at both 
qubit and continuous variable levels and, in Sec. |6] the violation of multipartite non locality 
tests, and show them to be within the reach of current technology. Radial modes will thereby 
turn out to be promising not only for information processing but also as probes of fundamental 
physics. 

2. The trap 

We shall consider the radial modes of n ions of mass m and charge ze in a linear Paul trap 
131. Let Xj and Pj be the position and momentum operators associated to the radial degree of 
freedom of the j -th ion, which is trapped in the radial direction with angular frequency Uj (see 
Fig [Til. In the following, the longitudinal trapping frequency lol will be the unit of frequency 
and will set the unit of length as well (equal to c? = ^ z^e^ / {A'K£[)muj\), where Eq is the 
dielectric constant); also, we shall set ?i = 1, so that all the quantities will be dimensionless. 
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Figure 1. Considered setup: the ions ai'e trapped along the longitudinal direction with the same 
trapping frequency cj^. Also, each ion j is trapped along the transverse, 'radial' direction (i.e., 
along the arrows in the drawing) with trapping frequency ujj . The operator Xj stands for the 
position of ion j along the radial direction: these radial oscillations are described by the modes 
{Xj , Pj } we shall focus on. 



The Coulomb interaction affects the local radial oscillation frequencies: for convenience, let 
us then define the 'effective' local radial frequencies 



\ 
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{uj} being the equilibrium positions of the ions O. Rescaling the canonical operators 
according to Xj = ^JmujjXj, pj = Pj/ y/mCjj , and grouping them in a vector of operators 

R = (ii,... , Xn,Pj , ■ ■ ■ , Pn) , allows onc to express the global Hamiltonian of the system 
in the harmonic approximatiorjj as 

H =^R^ {k®Cj)R, (1) 

where a; is a diagonal matrix: uj = diag {uji, . . . ,(!;„), while the potential matrix k has 
diagonal entries Kjj = ujj and off-diagonal entries Kjk = lo\(P /{yJujjUJk\ Uj — u^^) for 

Let us recall the canonical commutation relations \Rj ,Rk] = i^jk 1, where the 2n x 2n 
matrix fl has entries fij ^ = 6n,k-j — Snj-k for 1 < j,k < 2n, and that Gaussian states are 
defined as states with Gaussian characteristic functions: a Gaussian state g is thus uniquely 
determined by its "covariance matrix" (CM) cr, with entries given by 

a,k = ^Tr {{Rj,Rk)e\ - Tr {Rjq\Tt {RkQ\ , (2) 

and by the vector of first moments i?, with components Rj = Tr \Rj g] H |5l . 

The ground state of Hamiltonian is a Gaussian state with a block diagonal CM 
— {(Tx © (T~i)/2, where 

<T, =<I;l/2(<:,l/2^(I,l/2)-l/2^1/2 ^ 

and vanishing first moments. 

t Notice that, in the harmonic approximation (i.e., at second order in the ions' displacements) the couphng between 
radial modes and longitudinal ones vanishes. Fourth order couplings are ai'ound 10^"* times smaller in the considered 
experimental conditions and will thus be safely neglected. 
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Finally, let us remember that the evolution for the time t of an initial Gaussian state with 
CM a under a quadratic Hamiltonian H = W HR (where H is any symmetric matrix) is a 
Gaussian state with CM given by 

a-t = Sto-Sj , (3) 

for St = expiflHt). 



2.1. Dissipation 

In our study, we will take into account the decoherence of the radial modes in an environment 
of phonons with temperature T and 'loss rate' 7 (for simplicity assumed to be the same for 
every mode). Under such conditions, the evolution of the ions' state g at frequencies {Cjj} is 
described by the following master equation in interaction picture ||6l 



^ = ^ XI [^j(2a]£'aj " f^jO^^jQ - QO-jf^]) + {^3 + l)(2aj£'a] - a]a.jg- Qa]aj) 



(4) 



where the number of phonons in the radial mode j is given by A^j := 1/ ^exp (^^^ 

according to Bose law (ks being Boltzmann constant) and aj := {Xj + iPj)/\/2. In 
accordance with the experimental terminology, we will also refer to the quantity := Nj^ 
as the 'heating rate', essentially representing the rate at which thermal phonons are injected 
into mode j. 

Eq. IHpreserves the Gaussian character of the initial state and can results, for the CM cr 
of Gaussian states, into the following equation: 

d 

— cr = 7croo - icr 

(where cr 00 = cr' © cr', cr' = Aia.g{Ni + \, . . . , Nn + \)), with solution: 

e-^Vt + (1 - e-'^*)ffoo . 



3. Linear phononics and beyond 

In this section, we shall assume that the trapping frequencies {wj}, and thus {ojj}, can be 
controlled locally and changed suddenly {i.e., much faster than w~^). Such a local control 
may be achieved by building small, local radial electrodes, by adding local optical standing 
waveslH or in trap arrays EUD. Let us also note that a setting formally identical to the one 
we describe here could be reproduced shortly even for locally defined axial ('longitudinal') 
modes in segmented Paul traps ||9l. The changes in trapping frequency we will consider 
(typically ranging between 1 and 10 MHz) can be realised in about 10ns ifTol . and will be 
regarded as instantaneous over our timescales. Our aim is to show how, in principle, local 
control of the frequencies allows one to perform any arbitrary "linear optical" operation on 
the radial modes of the ions, that is any unitary operation under which, in the Heisenberg 
picture, the vector of operators R transforms linearly: R ^ SR. The matrix S has to be 
'symplectic', i.e. such that 

s'^ns = n 

X Currently, standing waves realize, at most, trapping potentials of about 1 MHz in experiments where the internal 
degrees of freedom are controlled. However, the manipulation of radial modes can tolerate much higher scattering 
rates than that of internal degrees of freedom, and would thus allow for much higher trapping frequencies. 
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(where is the anti-symmetric symplectic form, defined at page [3]), to preserve the 
canonical commutation relations. Notice that such a class of operators includes squeezing 
transformations H We should mention that the idea of squeezing the state of a single ion by 
controlling the trapping frequency was first put forward in the early 90's ifTTl . Also, single- 
mode squeezing can be achieved by coupling the motion of the ion with its internal degrees 
of freedom llT2HT3l . The scope of the present Section is however wider, since we set out to 
prove that any symplectic operation on any number of ions can be implemented by frequency 
manipulation. 

Let us first remark that any symplectic operation S* on a system of many canonical 
degrees of freedom ("modes") can be decomposed into a combination of generic single- 
mode symplectic transformations and two-mode rotations ("beam splitters", in the quantum 
optical terminology) lfT4l[T5l . This fact follows from the Euler decomposition of symplectic 
operations lIT?! and from the possibility of decomposing energy preserving operations into 
a network of beam-splitters between pairs of modes ITSl . It is therefore sufficient for 
us to establish the possibility of performing these two subclasses of operations (single- 
mode symplectic transformations and two-mode beam splitters) on our system of n ions by 
manipulating the local frequencies. In turn, again because of Euler decomposition, single- 
mode operations can always be reduced to combinations of squeezings and phase-shifts. 

3.1. Single mode operations 

In what follows we assume that the original effective frequencies of the ions are different 
but commensurate, as given by, say, Cjj = jcj, and that uj is large enough so that interaction 
between any two ions is suppressed when the relevant coherent manipulation sets off[j] Let 
us then consider the reaction of the system if the frequency of the j-th ion changes suddenly 
from Cbj to ajCjj, for some real Uj (clearly, such changes are completely governed by changes 
in the w/s). The Heisenberg equation of motion for Xj andpj can be immediately integrated 
in such a case, resulting into a symplectic transformation Sj{t) given by 

The first and last matrix of this decomposition correspond to 'squeezing' operations in the 
quantum optical terminology, whereas the second factor is known as a 'phase shift' {i.e., a 
rotation in the single-mode phase space). Combinations of squeezings and phase-shifts make 
up any possible single-mode symplectic operation: we thus need to show that such operations 
can be implemented individually on any ion of the system in a controllable manner (as pointed 
out above, this is sufficient because of the Euler decomposition of single mode operations 
Gl). 

§ In the literature, squeezing transformations liave a tendency to be excluded from the class of "linear" 
transformations, because their implementations require 'nonlinear' (third order) interactions. However, the effective 
evolution of the modes of interest ('signal' and 'idle', usually) is actually linear in the sense specified above (linear 
evolution of the vector R in Heisenberg picture), and we will thus generally refer to 'linear' transformations as 
including squeezing. 

II The linear scaling of the coupling is not a stringent requirement. Because the coupling between ion j and k falls 
off like \j — fcl""^, one can safely assume, say, a 'triangle' profile (alternately increasing and decreasing) for the 
couplings in long chains of ions. Note also that, for a linear scaling, the population of the levels of ion j due to the 
interaction with k roughly scales as 1/(1 + {uj/ijJi^)-^\j — fcp): a frequency step li ~ 20a; is already enough to 
make the effects of all the interactions essentially negligible. 
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3.1.1. Phase-shift To realise a phase-shift operation on the fc-th ion, without any squeezing, 
it is sufficient to change the frequencies of all the other ions in the same way, such that, in the 
notation defined above, a)^ ~ 1 and aj = a ^ \ for j ^ k {aj being the factor by which 
the frequency of ion j is multiplied by manipulating the local trapping potential). As apparent 
from Eq. (|5]l, after a time ta ~ 2-K/(uja) from the change one has Sj ~ I2 for j 7^ k: all 
the ions which undergo a change in trapping frequency are back to their initial state (let us 
recall that Cjj = jCj by assumption, so that after the central rotation of Eq. (|5]) reduces to 
the identity and the two opposite squeezing transformations just nullify each other). On the 
other hand, the transformation Sk will be equal to the central rotation (with no squeezing, as 
ak is kept equal to 1) of angle tp^. = 2'Kk/a: the oscillation of the fc-th ion will have acquired 
such a phase, analogous to the "optical phase" of light modes. If the frequencies are switched 
back to the original values after a time ta, the net effect of the evolution is then analogous to 
an 'optical' phase-shift of angle ipk on the ion fc. 

3.1.2. Squeezing In order to squeeze the state of ion fc, leaving all the other ions unaffected, 
one can conversely change only the pertinent frequency, so that au ^ \ and aj — 1 for 
j 7^ fc. Then, after a time period t = 2tt /uj, all the other ions will have returned to the initial 
state, while ion k will be squeezed and phase-shifted according to Eq. (|5]). The phase-shift can 
always be arbitrarily corrected by applying the strategy to obtain "pure" phase shifts described 
in Sec. 13.1.11 which conclusively shows that any linear operation on a single mode can be 
implemented by controlling the local trapping frequencies]^ 

The degree of squeezing achievable with such a strategy depends crucially on the phase- 
shift operation since, as shown by Eq. (|5]), the two squeezing operations act along orthogonal 
directions and are the inverse of each other. In the case ak = {j + h)/k for /i e IN, the two 
squeezings act effectively along the same phase space direction: in this instance the phase- 
shift can be balanced by a counter-rotation of tt/A in phase space and the final squeezing 
operation is a diagonal matrix given by diag(afc, «fe In principle, the value of ak is only 
limited by the stability of the system and the breakdown of the harmonic approximation 0, 
occurring when the squeezing is comparable to the ratio between the size of the wave packet 
and the distance between the atoms. In actual experiments, the considered setup would thus 
allow for Y^Ofc ^ 500. Notice that such values are by far out of reach in optical systems, 
where squeezings corresponding to ~ 8 were recently reported |fT6l . Also notice that, 
by placing the ions inside cavities, the squeezing of the massive degrees of freedom could be 
transferred to light, so that radial modes could act as an effective source of squeezing (and 
potentially even entanglement) for optical systems as well ITtI . 

3.1.3. Beam-splitters To extend our proof to any linear operation, on any number of ions, 
let us now turn to 'beam-splitting' operations between the radial modes of any two ions in the 
array. Such operations are achieved by bringing the two modes in question (hereafter labeled 
by j and fc) to the same frequency ui = coj = Cok, so that the Coulomb interaction between 
them is no longer suppressed. Switching to the interaction picture, one has the following 
coupling Hamiltonian between the two modes: 

Hi = KjkXj{t)xk{t) = K,k{a, e-*"^* + a] e*'^*)(afc c"''^* + 4 c^"*) , 

^ In fact, arbitrary adjustment of the phase allows for "pure" squeezing operations along the phase space directions 
X and p, and thus implies the possibility to construct any linear operation, because of the symplectic Euler 
decomposition previously discussed. 
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where the ladder operators are defined as Xj ~ {aj + a^) and Kj^ is an entry of the matrix 
K defined in Eq. ([T]i. If the frequency a) is sufficiently large the rotating wave approximation 
reliably applies and the rotating terms can be neglected to yielcS 

Hj = Kjkiajal + a]afc) . 

This Hamiltonian realises exactly the desired beam splitter-like evolution, resulting into a 
symplectic transformation which mixes Xj with x^. and pj with p^ (rotating such pairs 
equally, by the angle Kjkt). For instance, a '50 : 50' beam splitter is achieved after a time 
t ~ 7r/(4Kjfe). Since the interaction requires a change of the local frequencies it includes 
automatically in it a local operation, which may however be corrected before or after the 
'beam-splitting' procedure, as detailed above. 

Summing up, the combined arguments presented in Sees. I3.1T1 13. 1.21 and 13. 1.3] show 
that any "linear optical" operation, including squeezing operations and entangling operations 
acting on multiple different modes, can be implemented for radial modes of trapped ions by 
properly tuning the frequencies of the radial microtraps of the individual ions. 

Linear optical operations, complemented by displacements (see Sec. 13. 2t . correspond 
to all the unitary transformations that preserve the Gaussian character of the initial state. 
Therefore, all the developments based on Gaussian states in quantum optics, in particular 
concerning entanglement manipulation |4] and information protocols flS], can be carried 
over to radial modes of ion traps if local control is achieved. Arguably, the framework we 
have outlined for the implementation of linear optical operations on the local oscillations 
of distinct ions, to which one may refer as "linear phononics", could offer further practical 
advantages over its quantum optical counterparts, essentially related to the static nature of the 
quantum information processed (for instance, 'mode matching', which hampers substantially 
the implementation of linear optics, is not a problem here). 

3.2. Further manipulations and measurements 

Beside the linear operations treated so far, implemented through varying trapping potentials 
and Coulomb interactions, the motional degrees of freedom of trapped ions allow for other 
controlled manipulations, some of which take advantage of the coupling between motion and 
the internal degrees of freedom that can be engineered by applying standing- or traveling- 
wave pulses to the ions. In order to give a complete account of the possibilities offered by 
radial modes, let us here review such strategies, and briefly comment about their applicability 
to radial modes. As a general remark let us mention that, because of the tighter confinement 
they allow for, radial modes easily meet the Lamb-Dicke condition (depending on the width 
of the ground state's wavepacket), which means that the coupling with the internal degrees of 
freedom can be tailored to a high degree of accuracy (generally better than for longitudinal 
modes)13 Also, individual ions can be addressed in such manipulations as, at spacings of 
some micrometers and assuming pulses' waists of the order of 1 fim, less than 1% of the 
central laser power would be shined on neighbouring ions with respect to the central one. 

Displacement operations, which shift the operators Rj by real numbers, can be realized 
in several ways: by classical driving fields, by standing waves, or by shifting the radial 
equilibrium positions of the ions |fT9l . In particular, transverse driving fields could be 
applied to displace the radial modes, with no particular hurdles. In the following, the unitary 

This always holds at the trapping frequencies we shall consider in specific examples, going from 1 to few tens of 
MHz, for which stability around 10^^ are achievable. 

* The price to pay for such an accuracy is longer operation times, which could be ultimately reduced if stronger 
lasers became viable. 
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displacing the canonical operators of mode j by, respectively, Xj and pj will be denoted by 

bj{:xj,p^j). 

Notably, even non-Gaussian states can be engineered in this setup with relative ease 
(with respect to quantum optics, where non-Gaussian manipulations require higher order 
nonlinearities, usually extremely weak), either by entering the nonlinear regime of the 
Coulomb interactions or by coupling the internal degrees of freedom of the ions. The 
experimental realisation of a cross-Kerr coupling for the longitudinal oscillations of trapped 
ions has been recently reported 

The coupling to internal degrees of freedom also allows for Gaussian and non Gaussian 
measurements on individual ions. The measurement of quadrature operators, corresponding 
to homodyning, was proposed in 1201 1211 l22l . Quantum non-demolition measurement of 
local number states and parity could be measured as well, by applying the scheme suggested 
and realised for cavity QED in Ii231 , based on the dispersive coupling of the number of 
oscillations to two internal levels of the ion. Quite remarkably, such a scheme would allow 
one to measure the phonon number's parity on a single copy of the state and run of the 
apparatus. As we will see in Sec. |6] this possibility is consequential for the violation of 
Bell-like inequalities. Notice also that such non-Gaussian measurements pose a considerable 
technological challenge for light fields, where resolving photon numbers with high detector 
efficiency is still daunting despite recent progress |j24j,25J. On the other hand, the schemes 
recalled above are bound to be comparatively slow, requiring to wait for half a period of 
coherent interaction between motions and internal levels and then subsequent readout of 
the internal levels by, e.g., fluorescence. When internal degrees of freedom are involved, 
radial modes require longer times, roughly on the order of tens of microseconds, but achieve 
remarkable precision. 

Finally, concerning the preparation of initial states for coherent manipulations, let us 
remind that cooling of ion oscillations to their ground state can be achieved very efficiently 
by sideband-Raman pulses ll26ll27l (or more recent variations over such a strategy ll28ll29l ). 
Note, since in the described scenario the initial local potential is much larger than the 
interaction, simultaneous cooling of all the chain could be done by cooling each ion to its 
local ground state, which is to a good approximation the global ground state. The local single 
excitations on such ground states may be prepared by addressing the first ion with a proper 
sequence of blue and red sideband pulses, as detailed in lfT9l . 

In the remainder of the paper we shall demonstrate the potential of linear operations on 
radial modes in specific applications, with a particular focus on settings fully accessible to 
current experiments. 

4. Entanglement generation 

The set of operations described in Sec. |3] including local squeezing and two-mode beam 
splitters between distant modes, allows for the generation of entanglement between the radial 
modes of two ions in the chain. In fact, as is weU known, two single-mode states squeezed 
along orthogonal quadratures, entering a beam splitter, give rise to an entangled state for 
the outgoing modes (this is a standard procedure to create CV entanglement in quantum 
optics). In principle, for the system of trapped ions under examination, the local control 
of the transverse trapping frequencies allows for the creation of tailored entangled states, 
where the entanglement can be build up between any two ions in the chain. Moreover, 
and possibly even more intriguingly, the possibility of controlling a system of n parties (the 
ions) all constantly interacting with each other, paves the way for the creation of multipartite 
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Figure 2. Dynamics of the entanglement (in ebits of logartithmic negativity) between the radial 
modes of the ions, in a trap containing n = 4 ions and with longitudinal trapping frequency 
u)i^ = IMHz. At time i = the system starts off from the ground state for transverse 
trapping frequency uii = 5MHz and then evolves under a trapping frequency = 2MHz. 
The plots on the left show the logarithmic negativity shared by all the possible bipartitions of 
the four inodes (because of the spatial symmetry, one has only five distinct cases in the chain: 
for instance, one has E^^'^'^ = E^'^'^). On the right, the logarithmic negativity shared by 
pairs of inodes is plotted (again, because of syininetry, there occur only four distinct cases). 
The continuous (red) curves refer to a case with no dissipation, while the dashed (green) curves 
refer to a case with heating rate e = 2KHz (resulting from a coupling to the bath 7 = 10~*Hz 
and a temperature T = 294° K). 



entanglement between the distinct ions. 

In this section we show that unprecedented degrees of bipartite entanglement, as well 
as interesting and robust multipartite entangled states, can be obtained in settings accessible 
to current experiments, starting from a ground state of the system and requiring only global 
control of the trapping potential (where, i.e., the trapping frequencies {wj} are the same for 
all ions at any time, and are only changed simultaneously). Further, we will briefly discuss the 
possibilities a local control of the frequencies that would open up for the controlled dynamics 
of multipartite entanglement. 

The specific situation we shall address starts off from the ground state Qg of Hamiltonian 
H [see Eq. ([T]l] - with all frequencies being equal, i.e. Uj = uji for 1 < j < n - as the 
initial state. Next, the trapping frequencies are changed to the common value ujf, so that the 
state Qg will not be stationary anymore under the modified Hamiltonian. For large LOi, the 
initial state Qg contains very little entanglement but, if the frequencies are suddenly changed, 
entanglement builds up during the subsequent evolution (see |[30| for an analogous scheme 
in chains of nanomechanical oscillators). Entanglement between different subsystems will 
be quantified by the logarithmic negativity Ej^f = logj ||^||i, where ||^||i stands for the 
trace norm of the 'partially transposed' density matrix of the subsystem lISTl l5l. Different 
partitions will be denoted by superscripts referring to the ions: for instance, the logarithmic 
negativity between ions 1-2 and 3-4 will be denoted by Ej^'^'^. The intuition behind this 
entanglement generation method is that operating faster than the speed of sound in the system, 
i.e., essentially, faster than the inverse of the energy gaps, is analogous to operating locally. 
Thus, by changing the local potentials (in the same way for each ion) fast enough one 
generates local squeezing, which is then 'converted' into entanglement by the time-evolution 
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Figure 3. Dynamics of the entanglement (in ebits of logarithmic negativity) between the radial 
modes of the ions, in a trap containing n = 3 ions and with longitudinal trapping frequency 
u)L = IMHz. At time t = the systein starts off from the ground state for transverse 
trapping frequency uii = 20MHz and then evolves under a trapping frequency i^f = 2MHz. 
The plots on the left show the logarithmic negativity shared by all the possible bipartitions of 
the four modes (because of the spatial symmetry, one has only two distinct cases in the chain). 
On the right, the logarithmic negativity shared by pairs of inodes is plotted (again, because of 
symmetry, there occur only two distinct cases). The continuous (red) curves refer to a case with 
no dissipation, the dotted (blue) curves refer to a case with heating rate e = 200Hz (resulting 
from a coupling to the bath 7 = 10~^Hz and a temperature T = 294° K), while the dashed 
(green) curves refer to a case with heating rate e = 2KHz (resulting from a coupling to the 
bath 7 = 10~*Hz and a temperature T = 294°K). 



through the harmonic coupling between the ions, analogous to a set of beam splitters. 

The case portrayed in Fig. |2] represents four ions starting from the ground state for 
uji ~ 5MHz, which then evolves under the frequency w/ = 2MHz for (for a transverse 
trapping frequency ojl = IMHz). As can be seen, the global state gets entangled under 
any possible bipartition of the four modes, a situation which is referred to in the literature 
as "complete inseparability" ll32l . A more extensive inspection shows that this is a general 
property of the Hamiltonian at hand, even for larger number of ions. Of course, as shown 
in Fig. [3] complete inseparability also occurs for three ions: in Sec.|6]we shall see how the 
multipartite entanglement for three ions analysed in Fig. [3] could serve to violate Bell-like 
inequalities. Our plots also show that such an entanglement exhibits considerable resilience 
under realistic heating rates. In this respect, notice that the heating rate for the dashed (green) 
curved plotted in the graphs is e = 2KHz, and that heating rates as low as ~ lOHz have been 
recently observed in ion traps (even though for a single ion only ll33l[34l ): robust multipartite 
entanglement (between 5 and 10 cbits of logarithmic negativity) between the ions can thus be 
created. Such degrees of entanglement are inconceivable for multipartite systems in photonic 
systems where, furthermore, the manipulation of several modes tends to become awkward 
due to the increasing number of mode-matching conditions to be fulfilled. 

In the plots, we also report the entanglement of pairs of modes in the traps for four, 
three and two ions. For bipartite entanglement as well, as shown in Fig. H] the degree of 
robust entanglement achievable is by far beyond the maximal values obtained in quantum 
optical systems^ even considering rather modest frequency jumps. The case portrayed in 

fj To the best of our knowledge, the highest measured value for the logarithmic negativity in optical sy steins (infeiTed 
from state reconstruction) is ~ 1.6 ebits 1351 . A simple evaluation also shows that, exploiting the degree of 
squeezing reported in 1161 and assuming perfect beam-splitting operations, one could achieve at most Ej^ < 3 ebits. 
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Figure 4. Dynamics of the entanglement (in ebits of logarithmic negativity) between the radial 
modes of the ions, in a trap containing n = 2 ions and with longitudinal trapping frequency 
= IMHz. In plot (a), the system starts off from the ground state for transverse trapping 
frequency a;; = 20MHz and then evolves under a constant trapping frequency = 2MHz. 
In plot (b), the initial frequency is LOi = 20MHz and then the system evolves for 5 /^s at 
frequency ojj = 2MHz and for 25 at frequency LOi ; for iterations of the cycle are displayed. 
The continuous (red) curves refer to a case with no dissipation, the dotted (blue) curves refer 
account for a heating rate e = 200Hz, while the dashed (green) curves refer to a case with 
heating rate e = 2KHz (resulting from a couphng to the bath 7 = 10~^Hz and a temperature 
T = 294°K). 



Fig. [4(b)] Starts from the ground state for Vi — lOMHz. Afterwards, the frequency is switched 
to w/ = 2MHz for 5/iS and then back to uji for 25/is: such a cycle reaUzes a highly 
entangling operation because it is tailored such that the squeezings act always in the same 
phase-space directions and, in principle, can be iterated at will to obtain any desired degree 
of entanglement. In practice decoherence and dissipation will degrade such an entanglement 
after some iterations. However, as is apparent from the plot, stable and robust entanglement 
(around 10 ebits of logarithmic negativity) between two ions could be created under realistic 
conditions. 

Also, in the traps with more than two ions, the dynamics of the entanglement between 
pairs of modes displays a 'monogamous' behaviour 1361 : for four ions, when the entanglement 
between ions 1 and 2 fades, the entanglement between 1 and 3 starts raising and then is 
replaced by entanglement between 1 and 4, which in turn fades out before the revival of 
the quantum correlations between 1 and 2. The entanglement between 2 and 3 - strongly 
favoured by the fact that such ions are both nearest neighbours and at perfect resonance with 
each other (not true in general, as the Coulomb corrections on the trapping frequencies depend 
on the position in the trap) - is always positive but, significantly, its peaks and dips follows 
those of the entanglement between 1 and 4, as one would heuristically expect by virtue of 
monogamy. As reminded above, the reported values of uJi and oj / represent the 'bare' trapping 
frequencies, not taking into account the corrections due to Coulomb repulsion: therefore, the 
ions are not all at resonance during the evolution. Now, if one assumes that local control were 
available, one could correct for such a mismatch by applying different bare frequencies to ions 
in different positions along the trap's longitudinal axis. We applied such corrections in the four 
ions case considered above and found that, quite interestingly, one can spread correlations 
more evenly in this way and get closer to the ultimate bounds imposed by monogamy: after 
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lO^s (all the parameters being the same aside from the corrections) one ends up with a state 
where each radial mode is individually entangled with all the other radial modes (a situation 
which never arises without corrections, see Fig.|2l). 

The continuous variable entanglement between two radial modes of the ions could 
be swapped to light if cavities were added to the setup, with the potential to achieve 
unprecedented degrees of optical entanglement, as the parametric processes are definitely 
outperformed by the strategy outlined above based on the trapping frequency's control: this 
possibility will be the subject of a detailed investigation in a forthcoming paper iflTl . 

5. Propagation of quantum information 

The locally defined radial modes of the ions, interacting via the Hamiltonian ([B, are an 
example of harmonic chain, where the coupling between different oscillators can be, to a 
very good extent, controlled. The propagation of quantum information in such a setting 
has been proposed and analysed theoretically in 1371 . but not yet realised in practice since 
mechanical oscillators have yet to enter the quantum regime (due to the technical difficulties 
still encountered in controlling both the couplings and the decoherence and dissipation of 
chains of nano-oscillators 1381 ). Trapped ions could thus provide a very promising alternative 
to implement harmonic chains, realize quantum data buses |[39l and test the results predicted 
by theoretical studies. More generally, such a demonstration would be a further application 
of ion traps as quantum simulators, as already proposed in several past studies l40ll4ll . In 
this Section, we shall present an overview of the possibilities offered by the radial modes as 
harmonic chains for the propagation of quantum states, with quantitative investigations, in 
scenarios similar to those detailed in the preceding sections. 

5.1. Transmission of two-dimensional quantum states 

The transmission of quantum information encoded in finite dimensional quantum states (most 
often qubit states) through a chain of interacting quantum systems has drawn much attention in 
recent years l42l . as such systems are envisaged as possible quantum buses, linking different 
parts of future quantum processors. However, despite considerable efforts, finding systems 
where such a transmission can be realised with sufficiently small level of noise has proven to 
be very challenging. Therefore, good practical candidates for such tasks are still of interest. 
Proposals in this context have so far focused on the transmission of finite dimensional states 
in chains of interacting (effective or proper) 'spins'. Here, instead, we consider qubit states 
encoded in the bosonic Hilbert space of the radial modes of trapped ions and, guided by the 
results of the Section[3] we show that such modes are able to send qubit states through chains 
of trapped ions. 

As a preliminary remark, let us notice that, if the ratio between the radial and the 
longitudinal trapping frequencies is very large, LOj/ujL ^ 1 Vj', then the effect of the Coulomb 
interaction on the Hamiltonian governing the radial modes is neglig ibl^ and the Hamiltonian 
([T]) reduces to iJ ~ © <^^- That is to say, the radial oscillations of the individual ions 

are decoupled, with local normal frequencies approximately given by {uoj}. In this situation, 
the ground state of the system is just the ordinary, non-squeezed, vacuum: a Gaussian state 
with CM equal to the identity and vanishing first moments. For instance, for LOj = SOcjl 
Vj the relative discrepancy between the actual ground state and the vacuum is only lO^'^ in 
terms of covariances. Starting from such a ground state, we will consider a scenario where 

tfSee page[3]and notice that, quite remarkably, tlie ratios d/\uj — nfe| does not depend on the trapping frequencies 
nor on the masses of the ions but only on the total number of ions n (see (3) for details). 
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Figure 5. Leftmost histograms: probability amplitude of a single excitation (phonon) across 
the chain, if the evolution occur for frequencies uij = SOojl Vj after times (from top to 
bottom) f = 0, f = 375c<j~^ and t = 5600aj~\ Middle histo grams: probability amplitude 
of a single excitation (phonon) across the chain, if the evolution occur for frequencies 
ljJ\ = ojio = Sojjj and ojj = 50w^^ for j ^ 1, 10 after times (from top to bottom) 4 = 0, 
t = 2800a; and t = 5600<^^^ Rightmost histo grams (red): probability amplitude of a 
single excitation (phonon) across the chain, if the evolution occur for frequencies u)j = IOujl 
\/j after times (from top to bottom) t = 0,t = 625aj^^ and t = 640ljJ^^ . 



the sender has its qubit state encoded in the single excitation sector of the oscillations of the 
first ion {i.e., on the subspace spanned by the phonon states |0) and |1)) and wants to send 
the excitation through the chain to a receiver who owns the last ion. Notice that, as recalled 
in Sec. 13.21 a (non-Gaussian) single qubit state of this kind could be prepared, in practice, by 
coupling the motion of an ion to the internal degrees of freedom through a proper sequence of 
blue and red sideband pulses. 

Neglecting, for the time being, thermal noise (see below), one has that the state |0) 
can be sent perfectly through the chain. Therefore, the overall fidelity of the transmission 
depends solely on the fidelity with which |1) can be sent from the first to the last ion. The 
amplitude of such a transfer can be determined analytically by tracking the evolution of the 
initial operator a\ (creating the excitation in the first ion) in Heisenberg picture. Because the 
system is harmonic the evolution of a\ results in a linear combination of field operators, which 
can be determined with standard techniques. We could thus consider different situations. The 
histograms we plot show the probability amplitude P for finding a single excitation at different 
sites in the chain. The resulting transmission fidelity (averaged over the Haar measure of a 
single qubit space) is given by 1 — -i (1 — \fPY + | V^(l ^ V^), monotonically increasing 
with P. We will consider a chain of n = 10 ions and an initial ground state for trapping 
frequencies Uj = SOwl Vj. 
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Fig. |5(a)| shows the case where the trapping frequencies are left unchanged after the 
initial preparation. Because the interacting terms, while very small, are still present, the 
excitation propagates because the local mode of the first ion is not exactly a normal mode 
of the system. After an evolution time t = 375aj£^, 0.76 of the probability amplitude is 
transferred to the final state. However, the propagation is rather dispersive, as each mode is 
on resonance with all the others. 

In order to boost the probability of transmitting the excitation, sender and receiver can 
switch the local trapping frequencies of their respective ions from the initial value LOi ~ hQuJL 
down to a common value cj/: in this way they can realise the beam splitting operation at 
distance described in Sec. 13. 1.31 which, in principle, would allow for perfect swapping. This 
specific example thus also illustrates our previous general discussion and shows with which 
precision and over which operating times are linear operations actually possible over a chain 
of ions. As one can see from Fig. |5(b)[ where w/ = So^l has been assumed, the beam splitting 
operation is virtually perfect: all the probability amplitude is gradually transferred to the final 
ion, while the other ions are never involved in the process. This quality in the transfer comes at 
the expense of transfer time: over 10 ions the beam splitting operation takes roughly an order 
of magnitude more (t ~ 5600a;^^) than the imperfect transfer considered in Fig. |5(a)| This is 
simply due to the fact that, by keeping the middle ions on resonance, one takes advantage of 
the interaction between nearby ions, which are clearly stronger (as the interaction decays like 
the cube of the distance) and propagate the excitation through the chain. For given trapping 
frequencies, the time needed to achieve the beam splitter between radial modes at the two end 
of a chain of n ions scales very accurately as n?. Notice also, as a side remark, that since 
the Coulomb corrections to the local trapping frequencies are symmetric with respect to the 
longitudinal centre of the trap, the transfer between extremal ions is somewhat favoured in 
practice, as equal trapping potentials will result in equal effective trapping frequencies for 
such ions. 

Finally, in Fig. |5(c)[ we consider a 'compromise' between the two cases analysed above: 
all the frequencies are changed to the frequency lOcj^ obtaining, after a time t = 625a;£^, 
a transferred probability amplitude of 0.85. It is important to remark that the slower beam 
splitting operation, where the quantum information does not disperse through the chain at all, 
is not only more precise than the other two options examined, but also much more stable with 
respect to imperfections in the allowed interaction times. 

Let us now turn to inspect the effect of decoherence: clearly, thermal phonons are 
deleterious for the transfer of information encoded in single excitation sectors. The analytical 
estimates for the transfer times we just presented allow us to determine the restrictions on the 
heating rates that would allow such systems to transfer information in practical situations. 
Assuming a longitudinal trapping frequency of 1 MHz and a (less ambitious) chain of 5 
ions, one would realise a beam splitting operation (perfect transfer) between extremal ions 
in t ~ 1.4msec. In such a case, heating rates as low as e ~ 0.1 kHz would be needed for a 
coherent transfer to take place. For such heating rates, the action of the thermal phonons on 
the fidelity between initial and final state can be assumed to be linear resulting, after average 
over the Haar measure, into a mean fidelity ~ 1 — 3et/2 ~ 0.8. | Such heating rates are 
certainly demanding for arrays of traps, but not inconceivable, considering the heating rates 
achieved for single traps and the fact that current systems all operate at room temperature 

The effect of losses and thermal phonons for small e (so that et <C 1 over the interesting time-scales) can be easily 
estimated by letting the master equation Eq. (4) for A'^ ^ 1 act at first order on an otherwise perfectly transferred 
generic state a|0) + and by setting e = '~fN . The fideUty between initial and final states (corresponding to the 
overlap, for pure initial state) can be then determined and averaged over the Haar measure of a single qubit Hilbert 
space, to obtain the mean value ~ 1 — 3et/2, 
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Figure 6. Propagation of the entanglement (logarithmic negativity) between internal and 
motional degrees of freedom. The initial state is a Bell singlet between the internal degrees 
of freedom of the second and the third ions; the transverse trapping frequency is lOcJi,. The 
upper (lower) plot on the left column shows the logarithmic negativity between the internal 
(motional) degrees of freedom of ions 2 and 3. The upper (lower) plot on the right column 
shows the logarithmic negativity between the internal (motional) degrees of freedom of ions 
1 and 4. The dynamics is divided to three stages. In the first stage (for < * S: ^''^'ZjJ^^' 
the coupling between internal and motional degrees of freedom of ions 2 and 3 is on, and the 
entanglement is swapped the former to the latter. In the next stage (for 1.22^ < t < ^■'^Z^^' 
all the internal and motional degrees of freedom couplings are off, and the entanglement 
propagates through the chain via the motional degrees of freedom. Finally, for t > 6.2^, 
the the coupling between internal and motional degrees of freedom of ions 1 and 4 is on, and 
the final swapping takes place. 



(thus, in principle, two orders of magnitude could be gained on heating rates). 
5.2. The role of the internal degrees of freedom 

In the preceding section, while discussing the transfer of finite dimensional quantum 
information, we started off from a situation where a qubit state is encoded in the first two 
number states of the first ion. As already mentioned above, one of the most expedient ways 
to create such a finite dimensional state for the motional degrees of freedom is creating 
the desired state for the 'internal' degrees of freedom of the ion (embodied by its internal 
electronic levels) and then coupling them to the local radial modes to achieve the swapping. 
In order to provide the reader with a more comprehensive treatment and give him/her a flavour 
of the way the internal degrees of freedom enter in such dynamics, we briefly give account of 
a specific example where the internal degrees of freedom are involved throughout the whole 
transmission of quantum information. 

We shall consider a trap with four ions, with ojj = IQojl for 1 < j < 4. To keep the 
example down to earth, we will not assume any capability of local control nor any change 
during the evolution for the radial frequencies here. The only manipulation required will be 
the switching on and off of the coupling between the local phonons and the two levels of each 
ion where the initial qubit state is encoded. Denoting such levels \e)j and \g)j, and setting 
Cx.j \s){g\j + |5)(e|j, we shall assume the following interaction Hamiltonian iJj"* for 
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each ion j: 

which describes appropriately the coupling realised in experiments (see, e.g., {191 ). For 
simplicity, the interaction strengths gj will be all set to ojl (notice that, because the oscillators' 
frequency are set to IQlul, the interaction reduces, up to a very good approximation, to a 
rotating wave one). As initial state, at t = 0, we assume a Bell state between the internal 
degrees of freedom of ions 2 and 3, and the ground state for the remainder f the system: 

l5)i «) (|<?)2 ® |e>3 + |e>2 «) |.g>3) «> {9)4 «> Qg , 

where Qg is the ground state of the radial modes of the ions. As a signature for the transmission 
of quantum information, we will consider the evolution of the entanglement between the 
internal and motional degrees of freedom of ions 2 — 3 (initially entangled) and ions 1 — 4 
(which gets entangled through the process), which is reported in Fig. |6]in tems of logarithmic 
negativity (the logarithmic negativity between internal degrees of freedom of ions j and k is 
denoted E-j^). The dynamics is then split into three stage: from time t = to i = 1.2 • 27t/ujl 
the internal and motional degrees of freedom of ions 2 and 3 are coupled and the entanglement 
is dynamically swapped from the former to the latter. Next, from t = 1.2-27r/(jJito 
t = 6.2 • 2tt/u!i^, the coupling is switched off and the entanglement propagates through 
the chain via the motional degrees of freedom, and the radial modes of ions 1 and 4 gets 
entangled. Finally, from t = 6.2-27r/ti>L to t = 10-27r/ti>L, the entanglementis swapped to the 
internal degrees of freedom of ions 1 and 4. One can see that, though the quantum information 
partially disperses under such conditions (essentially due to the stringent restrictions we put, 
in this instance, on the control of the local frequencies), the process is capable of transferring 
the entanglement from the internal degrees of freedom of ions 2 and 3 to those of ions 1 and 
4 through the radial motional degrees of freedom. 

5.3. Propagation of continuous variable states and entanglement 

Clearly, quantum information (and entanglement) can be propagated through the chains 
under examination also at a continuous variable level, when populating the whole infinite 
dimensional Hilbert space. The study of information and entanglement propagation over 
harmonic chains has been proposed and discussed in detail in ifJTl |39l . Here, we shall 
investigate cases of such propagations along the chains of ions we are considering, addressing 
specific instances which can be realised in the laboratory. Our aim is pointing out at what 
degree and under what conditions can the theoretical schemes based on harmonic chains be 
implemented on radial modes of trapped ions. 

As argued in JSTj, the capacity of transmitting quantum states between distant parties is 
closely related to the capacity of transmitting, or "swapping", entanglement between them. In 
fact, the latter also critically requires the transfer to happen coherently through the whole 
process. In view of this fact we will limit ourselves to consider the propagation of CV 
entanglement (rather than the fidelity between sent and received quantum states). Also, in 
order to focus on a feasible scenario, we shall consider a common initial condition: the 
chain of n ions starts off from the (completely separable) ground state for Uj ~ Vj and, 
then, entanglement between the first and the second ion is created, as described in Sec. |4] 
by switching their frequencies to the same frequency ajg. From now on the first ion is left 
off-resonance at frequency uj,, (and thus effectively isolated from the chain), and the aim is 
transmitting its initial entanglement with the second ion through the chain over to the final 
n— th ion. Two options will be examined: 
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Figure 7. Transfer of the entanglement with the first ion of a chain (measured by ebits of 
logarithmic negativity i?^) from the second ion to the last one, in a chain of four ions. From 
top to bottom, the plots show the entanglement between ions 1 and 2 (E^), 1 and 3 (E^) 
and 1 and 4 (Ej^). On the left side (a), the entanglement is transferred by a relay of beam 
splitters obtained by switching the radial trapping frequencies of the ions 2 and 3 to \/98.8lui^ 
for a time 46t<j£ ^ while keeping the other ions off-resonant and then, finally, by switching the 
trapping frequencies of ions 3 and 4 to VSS.Sa;^ (with the other ions off-resonant). On the 
right side (b), the entanglement is transferred directly by a beam-splitter between ions 2 and 
4, obtained by setting their frequencies to \/98.Sujl while keeping the others off resonance. 
In both cases, the state at time i = is the ground state for all trapping frequencies equal to 
\/98.8ujl (such a value takes into account the corrections due to Coulomb repulsions) and, 
duiing the initial 2 /is, the entanglement between ions 1 and 2 is built up by bringing their 
'bare' trapping frequencies down to 2 fis. 



i) the entanglement is swapped by a chain of beam splitters between neighbouring ions 
(achieved by setting such neighbouring ions at the original frequency while putting the 
others off-resonance); 

ii) the entanglement is swapped directly by a beam splitter between the ions at the far ends 
of the chain. 

Notice that in case i) the beam-splitters operate at the original frequency luq so that no local 
squeezing takes place, as evident from Eq. (|5]l. Local squeezing would affect the entanglement 
between the modes but also, crucially, alter the original state to be transmitted and is thus 
undesirable in the present context. 

The two strategies i) and ii) are compared in Fig. O for a chain of four ions. Such a small 
system already permits one to highlight all the essential features of the two strategies. As 
apparent, a complete dynamical swapping of the entanglement can be achieved in both cases, 
with remarkable accuracy. A chain of beam splitters allows for a faster transfer, as one should 
expect since the coupling Kjf^ between distant ions is inversely proportional to their cubed 
distance: in general, the time needed to send information across n ions by a 'relay' of beam 
splitters roughly scales like n, whereas the time needed to achieve a beam-splitting operation 
between first and last ion scales like n^. A 'relay', taking less time, is thus less sensitive to 
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decoherence and dissipation. On the other hand, a single beam-splitting operation between 
distant ions only requires a single adjustment of the trapping frequencies (after the initial 
entanglement is created), whereas a chain of beam-splitters requires ?? — 2 manipulations, 
which would be possibly difficult to master in practice and could involve more errors and 
imperfections. Therefore, we have investigated the effect of imperfections in frequencies and 
evolution times on the transfer for both cases: the series of beam-splitters turns out to be more 
robust. In the instance depicted in Fig. |7] the transfer by a series of beam-splitters is left 
virtually unaffected by uncertainties of QAlol on frequencies and O.loj^^ on operating times 
{ujl is the longitudinal trapping frequency) whereas, on average, the entanglement transferred 
by a single beam splitter is reduced to to ~ 1.2 ebits of logarithmic negativity (out of ~ 2.2 
ebits initially present between ions 1 and 2) by the same imperfections. Summing up, strategy 
i) proved to be more reliable and faster (and thus less subject to decoherence and dissipation) 
and is hence to be preferred for the transmission of CV quantum information over chains of 
ions. Considering once again the example of Fig. |7] for a longitudinal frequency ujl ~ \ 
MHz, an effective coherent transfer would take about 100 and thus require heating rates 
around 1 KHz to be carried out effectively. 

6. Nonlocality tests 

The presence of strong Gaussian multipartite entanglement highlighted in Sec. |4] can be 
experimentally demonstrated and put to use in testing quantum nonlocality. Central to this 
endeavour is the capability of performing non-Gaussian measurements on the motional state 
of the ions, pointed out in Sec. 13. 21 In particular, we already recalled that phonon-numer parity 
measures are possible on single copies of the system, and so are displacement operations. In 
this section, we will explore the possibility of violating multipartite Bell-like inequalities 
(the so-called "Klyshko" inequalities ll43l ) by measuring "displaced parity" observables, as 
proposed for generic CV systems by Banaszek and Wodkiewicz B4l . Although bound to be 
subject to the locality loophole (considering that the distances between the ions are typically 
of the order of Ifim), such an experiment would be a remarkable test of quantum non-locality 
with massive particles, which is still lacking. 

To fix ideas and address a situation within the reach of current experiments, we shall 
study the test on the three-mode Gaussian state whose entanglement is described in Fig. [3] 
setting the evolution time to t = 5/iS and the heating rate to e = 200Hz (resulting from 
a coupling to the bath 7 = lO^^Hz and a temperature T = 294°K), whose CM will 
be denoted by 0-3. The family of (non-Gaussian) local, bounded, dichotomic observables 
for the displaced parity test is given by Ilj{xj,pj) = Dj{xj,pj)'^ {—ly''^ Dj{xj,pj), 
where Dj and fij are respectively the displacement and number of phonons operators of 
ion j. The three observers, pertaining to the three ions, randomly apply two different 
displacements [Dj{xj,pj) and Dj{xj,p'j)] on their ions and then measure parity locally. 
Such a measurement is clearly non-Gaussian, and allows one to violate Bell inequalities with 
Gaussian states. The expectation value of the operator H(i?) = Hi(a;i,pi) (8) Il2{x2,P2) <8) 
H3(a;3,p3) is simply proportional to the Wigner function W{R) of the composite system 
evaluated in the point i? = {xi, X2, X3,pi,p2,P3y : (H(i?)) = {2/Tr)^W{R) 145]. For a 
three-mode Gaussian state with covariance matrix cr^ one has 

WiR) = c-5flV^'^/(7rVDet 0-3) . 

The Bell-Klyshko inequality finally reads: 

g 

B3 = -^\W{xi,X2,X3,pi,P2,P3) + W{xi,x'2,X3,Pl,p'2,P3) 
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Figure 8. Function for Pi = P2 = P3 = p'l = P2 ~ P3 ~ 0' ^'1 ~ ^'3 ~ 6nm, 
X2 = — 4nm, — 7nm < x\ = < 3nm and — 3nm < X2 < 7nm for the three- 
mode Gaussian state obtained as detailed in the caption of Fig. [5] after a time t = 5fj,s. On 
the left hand-side no dissipation is accounted for, while on the right-hand side a heating rate 
e ~ 200Hz is introduced. Dimensions were reintroduced assuming uj^ = IMHz and Ca 
ions. The plane B3 = 2, above which violation occurs, is the bottom of the plots; in the 
displayed region, the functions reach maxima of ~ 2.45 for the pure state case and of ~ 2.28 
with dissipation. 

+ W{x[,X2,X3,p[,P2,P3) -W{x[,X2,x'ri,p[,P2,P3)\ < 2. (6) 

Quantum mechanics allows for 2 < i?3 < 4. 

Fig-Eshow a region in the space of displacements where the violation of the inequality 
is substantial and remarkably stable. It is apparent that, for such choices of displacements, 
the tolerable error on the displacement operation needed to maintain the violation is around 
1.5nm, which is within reach of current experimental capabilities. Also, as shown in the plot 
on the right, heating rates of 200Hz are still compatible with a violation of the inequalities. 

This preliminary study reveals very promising perspectives concerning the violation of 
Bell inequalities with massive degrees of freedom. Even if subject to a locality loophole, this 
endeavour would still stand out as a major, not yet probed, testing ground for fundamental 
quantum mechanics B6ll . and epitomises the considerable promise radial modes hold for 
quantum information and fundamental investigations alike. 

7. Summary and outlook 

We demonstrated how the local control of the trapping frequencies would allow one to 
reproduce any linear optical manipulation on radial modes of trapped ions. Drawing from 
previous studies on similar settings, we pointed out that phonon detection and homodyne 
detection as well as the implementation of non-Gaussian operations is possible in this setting. 
Next, we emphasized that, even restricting to global control, such manipulations enjoy a high 
efficiency in entanglement generation and a considerable resilience in the face of currently 
achievable dissipation rates. We then made manifest that radial modes could be used for the 
transmission of quantum information, stored in finite dimensional subspaces of the bosonic 
Hilbert space as well as in the full continuous variable domain. Finally we showed that, 
through (achievable) measurements of the displaced parity operator, Bell-like inequalities can 
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be violated under realistic decoherence and dissipation rates. As a rule of thumb, applying to 
systems with n ~ 6 ions or less and longitudinal trapping frequencies a;^ ~ 1 MHz, our study 
indicates that heating rates around 1 KHz are sufficiently low for coherent manipulations and 
robust entanglement generation in the continuous variable regime, whereas non-locality tests 
and transmission of finite-dimensional quantum dimensional information (stored in single 
phonons) are more delicate, requiring heating rates around 100 Hz to be carried out efficiently. 

In the light of the above, the experimental pursuit of the presented programme holds 
considerable promise, concerning both technological developments, such as the storage and 
manipulation of quantum information or the efficient generation of entanglement, and tests 
of fundamental physical aspects, as in the nonlocality test for massive degrees of freedom 
here discussed. As a first step to push this analysis further, a study is currently under way 
to examine the possibility to swap the remarkable entanglement produced in the traps to 
light modes, to ultimately beat parametric processes in the generation of quantum optical 
CV entanglement ini. 
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